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ABSTRACT 

We report our 110 ks Chandra observations of the supernova remnant (SNR) 0104-72.3 in the Smah 
Magellanic Cloud (SMC). The X-ray morphology shows two prominent lobes along the northwest- 
southeast direction and a soft faint arc in the east. Previous low resolution X-ray images attributed 
the unresolved emission from the southeastern lobe to a Be/X-ray star. Our high resolution Chandra 
data clearly shows that this emission is diffuse, shock-heated plasma, with negligible X-ray emission 
from the Be star. The eastern arc is positionally coincident with a filament seen in optical and infrared 
observations. Its X-ray spectrum is well fit by plasma of normal SMC abundances, suggesting that it is 
from shocked ambient gas. The X-ray spectra of the lobes show overabundant Fe, which is interpreted 
as emission from the reverse-shocked Fe-rich ejecta. The overall spectral characteristics of the lobes 
and the arc are similar to those of Type la SNRs, and we propose that SNR 0104-72.3 is the first 
case for a robust candidate Type la SNR in the SMC. On the other hand, the remnant appears to be 
interacting with dense clouds toward the east and to be associated with a nearby star-forming region. 
These features are unusual for a standard Type la SNR. Our results suggest an intriguing possibility 
that the progenitor of SNR 0104-72.3 might have been a white dwarf of a relatively young population. 
Subject headings: ISM: individual objects (SNR 0104-72.3) — ISM: supernova remnants — X-rays: 
ISM 



1. INTRODUCTION 

Among about a dozen supernova remnants (SNRs) de- 
tected in the Small Magellanic Cloud (SMC), only three 
X-ray brightest SNRs (0102-72.3, 0103-72.6, and 0049- 
73.6) have been relatively well studied. Because all these 
three SNRs most likely originate from a core-collapse 
SN, revealing the type of other SNRs is of great impor- 
tance to study the star-forming history, SN rates, and 
chemical evolution of the SMC. SNR 0104-72.3 is the 
fourth brightest X-ray SNR in the SMC, whose origin 
has been controversial. Its optical emission was sug- 
gested to be Balmer-dominated (Mathe wson et al.l ll984). 
indicating Ty pe la origin. On the oth er hand, with the 
ROSAT data, iHughes (fTool identified an un- 

resolved X-ray source in the southern part of the SNR 
as a candidate Be/X-ray star which appeared to be co- 
spatial with the SNR. Based on this plus other evidence 
for a Population I environment, they suggested a core- 
collapse orig:in for SNR 0104-72.3. With the XMM- 
Newton data, Ivan der Hey den et all ()2004 f) found evi- 
dence for enhanced Fe L line emission from the remnant, 
with which they reinstated Type la origin. Recently, 
with Akari IRC observations, Koo et al. (2007) discov- 
ered bright infrared shells which are positionally coinci- 
dent with th e optical Ha filaments surroundings; the X- 
ray emission (|Hughes fc Smithlll994f ). IR color distribu- 
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tions indicated shock interaction with ambient molecular 
clouds, which may support a core-collapse explosion of a 
massive progenitor star for the origin of SNR 0104-72.3. 

We report here on the results from our observations 
of SNR 0104-72.3 with the Chandra X-Ray observatory. 
We note that the remnant has been serendipitously de- 
tected a number of times by Chandra during the ACIS 
calibration observations of SNR 0102-7219. However, 
SNR 0104-72.3 was detected mostly in FI CCDs at large 
off-axis angles. Due to the low sensitivity and the poor 
spatial resolution, these archival data were not suitable 
for our analysis and we only report results based on our 
new observations. 

2. OBSERVATION & RESULTS 

We observed SNR 0104-72.3 with the Advanced CCD 
Imaging Spectrometer (ACIS) on board Chandra on 2008 
Jan 27-31 during A09. We used ACIS-S3 for the effec- 
tive detection of the soft X-ray emission of the SNR. A 
total of 110 ks of exposure, with two separate obsIDs (54 
ks for obsID 9100, 56 ks for obsID 9810) was obtained 
in the VFAINT mode. The level 1 event files were re- 
processed to create new level 2 event files. We applied 
parameters of the standard Chandra pipeline process, 
except that we applied the VFAINT mode background 
cleaning and turned off the pixel randomization. CIAO 
4.2 and CALDB 4.2.2 were used for all the reprocess- 
ing and analysis. We examined the overall background 
light curve for periods of high background. No significant 
background flare was seen and we did not apply any light 
curve filtering. 

3. X-RAY IMAGES 

Figure [TJa) is the broad band (0.6-3.0 keV) Chan- 
dra ACIS image of SNR 0104-72.3. Our deep ACIS 
image reveals details of the X-ray structures that 
have not been available from the previous ROSAT 
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Fig. 1. — (a) Chandra image (0.6-3.0 keV) of SNR 0104-72.3 smoothed with Gaussian beam of 2.5^^ (inset) zoom-up of the bright 
southern lobe, where previous observations suggested as a point source, smoothed with V beam. The location of the candidate Be star is 
mark ed with a white plus sign, (b) RGB composite image of the Chandra data overlaid with white contours of the 843 MHz radio map 
(^Bock et al.lfT99 9). Red is 0.6-0.95 keV, green is 0.95-1.12 keV, and blue is 1.12-3 keV bands. 



and XMM-Newton observat ions (jHughes fc SmithI 119941 : 
Ivan der Hevden et alJ[2QQ3 ). The X-ray morphology of 
the remnant is dominated by bright lobes extending 
about 90^' along the northwest-southeast direction. Faint 
emission surrounding the lobes is also clearly seen. Most 
notably, there is a faint arc-like structure on the east- 
ern side of the lobes. We refer to the region including 
two lobes as the "bar" region, and the faint arc emission 
in the east as the "arc" region (see inset image of Fig- 
ure [2l^a)). Figure [TJb) shows a RGB composite image of 
SNR 0104-72.3. The bar and the arc regions show dif- 
ferent colors in the RGB image. The arc is dominated 
by red while the bar shows a complex mixture of three 
colors. A color variation is also noticed within the 
bar, where the bright southern and northern regions of 
the bar are more greenish (i.e., 0.95-1.12 keV) and/or 
blueish (i.e., 1.12-3 keV) than the regions in the middle. 
There is thus a trend of higher photon energies toward 
the ends of the bar than in the middle. 



4. SPECTRAL ANALYSIS 

The high-resolution ACIS imaging observations allow 
us to perform a spatially resolved spectral analysis of 
SNR 0104-72.3. We extracted X-ray spectra from sev- 
eral characteristic regions, based on the color variation 
in the RGB composite image. For a background sub- 
traction, we used an average spectrum from four circular 
regions (^ 30'' in radius) around the SNR. For spectral 
modeling, unbinned spectra we re fit using the Churazov 
weighting (Churazov et al."1996). 

We first extract spectra from the arc and the bar. The 
background-subtracted spectra are shown in Figure [2l 
The observed spectra are soft, with most of the emis- 
sion seen below 3 keV. The peak energy of the arc is 
lower than that of the bar, which explains the differ- 
ent color between the arc and the bar in Figure [IJb). 
We fit the spectra with a single noneq uilibrium ioniza- 
tion (NEI) plane-parallel shock model (jBorkowski et al.l 
120011. vDshoc k in Xspec vl2) that is based on ATOMDB 
(|Smith et al.l l2001l) . We use an augmented version of 
this atomic database to include inner-shell processes and 



updates of the Fe L-shell lines Q We modeled the fore- 
ground absorption by the Galaxy and the SMC sepa- 
rately. For the absorption by the Galaxy, We used a fixed 
column density of A^H,Gai = 2.2 x 10^^ cm~^ which is the 
H I column density of the G alaxy toward SNR 0104- 
72.3 (I Pickev fc Lockmanlll990f ). The foreground column 
density in the SMC (A^h,smc) is fit as an additional ab- 
sorption component assuming the SMC interstellar abun- 
dances (Russell fc Dopita 1992). We first try to fit 
the spectra with elemental abundances of the emitting 
gasfixed at those of the SMC values (Russell fc Dopi taj 
1 19921 ). An acceptable fit is obtained for the spectrum of 
the arc region. The best fit results give hydrogen col- 
umn density of A^h,smc ~ 1 x 10^^ cm~^ and an electron 
temperature of 0.54 keV, although the ionization time 
scale is poorly constrained. This is the region where the 
remnant is thought to interact with dense ambient gas 
(see § [5]), and the X-ray spectral fit result is consistent 
with such an interpretation, i.e., the X-ray emission of 
the arc is predominantly from shocked gas with SMC 
abundances. On the other hand, the spectrum of the bar 
region cannot be adequately fit with models using SMC 
abundances (reduced chi-squared ^ 3); the prominent 
peak near 1 keV, in particular, could not be reproduced. 
The strong lines in this energy range are those of Ne K 
and Fe L shell lines. Freeing abundances of Ne and Fe 
provides an acceptable fit. The best fit model has over- 
abundant Fe (by a factor of 10 increase above the SMC 
value) and negligible amount of Ne. We note that quan- 
titative interpretations of the Fe abundance are difficult 
due to several modeli ng; issues including; in complete Fe 
L line modeling (e.g., "Badenes et al.|[2Q06[ ). and we in- 
terpret this number only in general sense of the Fe-rich 
ejecta. Our results of the Fe overabundanc e is co nsistent 
with the results of .van der Hevden et afl (|200l . The 
H I column density o f the SMC toward SNR 104-72.3 
is - 6 X 10^^ cm-2 C Stanimirovic et al.l 11999^ While 
the estimated A^h,smc ~ 9 x 10^-^ cm~^ is slightly larger 
than the H I column density, we note that HI column 

^ The augmented APEC atomic data have been provided by K. 
Borkowski. The relevant discussion on these detailed plasma model 
issues has been presented in Badenes et al. ( 2006). 
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Fig. 2. — (a) Chandra spectra of the bar and the arc regions. The inset image shows the extraction regions of the spectra. The sohd hnes 
are best fit models in Table [T] The spectra are binned for display purpose only, (b) Chandra spectra of subregions of the bar. The inset 
image shows the extraction regions of the spectra. The spectra of bar02, bar03, bar04 are scaled by lO"-*^, 10~^ and 10~^, respectively for 
a display purpose. The solid lines are best fit models in Table [T] scaled accordingly. 



density values typically underestimate the total hydrogen 
column density. 

The RGB composite image in Fig. [IJb) shows varying 
colors within the bar region, indicating spectral varia- 
tions although some of them may be simply due to low 
photon statistics. Figure Efb) shows extracted spectra 
from 4 smaller regions inside the bar, and spectra of 
some regions (e.g., bar04) show slightly different spec- 
tra than others. We fit spectra from each subregion with 
the vpshock model. Because of the poor photon statis- 
tics of individual regions, we could not constrain model 
parameters with both Ne and Fe abundances varied in 
the fits. Thus, as the previous fit results of the bar sug- 
gested negligible contribution of Ne, we fix Ne at the 
SMC abundance (fixing it at the best fit parameter of 
the bar region has only minor effect on the results). We 
further fix the absorbing column density at the best fit 
parameter of the bar region (A/'h,smc = 9-1 x lO^^cm"^). 
This approach provides an acceptable fit with reasonable 
parameter constraints. We find that the spectra from 
the central regions tend to have lower temperature and 
lower Fe abundance. The origin of this spectral variation 
is unclear with the current data. We speculate that the 
thermal conditions and/or metal abundances may be in- 
homogeneous along the bar region, possibly due to the 
mixture of the eject a and ambient medium and/or in- 
trinsically non-uniform distribution of the ejecta. And 
the spectral fits suggest that the central regions may be 
less dominated by ejecta material. All the fit results are 
tabulated in Table [H 



5. DISCUSSION 



5.1. 



The SN type of SNR 0104-72.3 

Based on t he ROSAT HRI observation, 
iHughes fc SmithI (|1994[ ) attributed the unresolved 
X-ray emission feature at a20oo, ^2000 = 01^ 06"^ 18®, 
-72° 06' 00'' to emission associated with a candidate Be 



star, interpreting it as a Be/X-ray star. The arcsecond 
resolution and good photon statistics of our deep 
Chandra data reveal that this X-ray emission feature is 
extended (Figure [1]), and no point-like source is found 
within a few arcsecond from the position of the Be 
star candidate (on-axis astrometric uncertainties of the 
ACIS is ~ 0.6"). We examined the RGB composite 
image and various narrow band images of different 
energy ranges and their ratio images, but could not 
find any point-like feature that is positionally associated 
with the Be candidate. The overall X-ray spectrum 
of this region is also similar to X-ray spectra in other 
regions. Thus we conclude that the X-ray emission 
in this region originates from the shocked SNR gas, 
and not from the poin t sour ce. Although the Be star 
that IHug hes fc SmithI (|1994D identified is unlikely to 
be a significant X-ray emitter, the possibility that 
SNR 0104-72.3 and the Be star belong to the same OB 
association, as they also suggested, remains valid, but 
would then argue against a typical Type la remnant. 
This will be further discussed in § 15.21 

The overabundance of Fe in the bar region suggests 
that SNR 0104-72.3 is likely the remnant of Type la SN. 
Indeed, the spectra of the bar regions are quite similar 
to ejecta s pectra found in the well-known Type la SNR 
DEM L71 (IHughes et al.ll2003D in the Large Magehanic 
Cloud (LMC) an d several other relative ly old Type la 
LMC SNRs (e.g., "Borkowski et ajjUool and references 
therein). DEM L71's central X-ray emission, which is 
well described by shocked Fe-rich ejecta, is similar to that 
of SNR 0104-72.3. DEM L71 also shows faint soft emis- 
sion of normal LMC abundances surrounding the central 
emission, similar to the arc in SNR 0104-72.3. Therefore, 
based on the overabundant Fe from our spectral fit and 
the spectral similarity of the bar and the arc region to 
other Type la SNRs, we propose that SNR 0104-72.3 is 
the first solid candidate for Type la SNR in the SMC. 



4 



Lee et al. 



TABLE 1 
Fit Results 



Region A^h,smc 


kT 


Ne 


Fe 


logr 


norm 


reduced 


[1021 cm-2] 


[keV] 








[10-4] 





Arc 9.9+5:9 

Bar 9.ll}j 

BarOl 9.1 

Bar02 9.1 

Bar03 9.1 

Bar04 9.1 



Note. — The errors are shown in 90% confidence range. Unless indicated with uncertainty 
values derived from the fits, the metal abundances are fixed at those of the SMC values from 
IRussell Domtal (1199^ ). They are He = 0.05, C = 0.13, N = 0.20, O = 0.10, Ne = 0.20, Mg = 
0.12, Si = 0.18, S = 0.15, Ar = 0.08, Ca = 0.20, Fe = 0.13, Ni = 0.20. All values are relative to 
the solar ones. The A^h,SMC values for the individual bar regions (BarOl, Bar02, etc) are fixed at 
the value obtained for the entire bar. 
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Fig. 3. — (a) Infrared image of SNR 0104-72.3 (IRAC 8/im of Spitzer archival data) with the contours of Chandra broadband image in 
Figure^a). (b) Ho; image from Hughes &; Smith (1994). The lowest level contour line from (a) is shown, (c) Ho; image around SNR 0104- 
72.3 from the Southern Ho; Sky Survey Atlas ( Gaustad et al. 2001 ). The solid circles are locations of known SNRs. The dashed circle is 
the approximate extent of the supperbubble DEM S124. 



Another element expected to b e a sig nature of Type la 
SNR is Si (e.g., Badenes et dl}\ 2007h . Allowing the Si 
abundance to vary in the fit yields a slightly larger value, 
but still comparable to the SMC abundance within the 
fit uncertainty. The Type la origin is supported by 
non-detection of a pulsar and/or pulsar wind nebula, al- 
though they could be too faint to be detect ed as in other 
core-collapse SNRs in the SMC fe.g.. IPark et al.ll2QQ3[ ). 

5.2. Type la SNR in Unusual Environment? 

Our deep Chandra observations find that the X-ray 
emission from the bar region is prominently from Fe- 
enriched SN ejecta of a Type la SNR. Based on limited 
sample, the X-ray morphologies of young Type la SNRs 
have been claimed to be statistically more symmetric 
than core-collapse SNRs (|Lopez et al.ll2QQ9l ). Although 
a quantitative analysis of the asymmetry in the mor- 
phology of SNR 0104-72.3 should be performed (which 
is beyond the scope of this work), the highly elongated 
morphology of SNR 0104-72.3 ap p ears t o be inconsistent 
with the results bv iLopez et al.l (|2009[ ). Such an elon- 
gated Fe-rich ejecta feature is similar to that found in 
the Galactic SNR W49B. Although the type is uncertain 
for W49B (e.g., Miceh et al. 2006), the highly elongated 
morphology of ejecta may suggest an asymmetric explo- 
sion. Alternatively, given the possibility that the rem- 
nant is interacting with dense material toward the east 
(see below), the atypical morphology could be due to 
highly inhomogeneous ambient environment of the rem- 



nant. 

The environment around SNR 0104-72.3 is also un- 
usual for a Type la SNR. The relatively high absorbing 
column density toward SNR 0104-72.3 indicates that the 
remnant could be located at dense environment. The 
eastern X-ray arc is positionally coincident with an op- 
tical and infrared shell (Figure [3]). The optical emission 
shows dominant Ha line with weak lines of other ions 
such as [O III] and [S II], and Payne et al. (2007) found 
that their line r atios are consiste nt with radiative shocks 
in other SNRs. iKoo et all (|2007[ ) reported infrared emis- 
sion associated with SNR 0104-72.3 from their Akari ob- 
servations. The observed infrared colors were consistent 
with t hose from shocked molecular clouds, and Koo et al.l 
(|2007[ ) suggested that SNR 0104-72.3 is probably inter- 
acting with molecular clouds. In this regard, the radio 
continuum emission of SNR 0104-72.3 could be mostly 
from the shocked ambient gas. Although the radio data 
quality is poor. Figure [TJb) shows that the extent of the 
radio continuum covers the eastern arc region, and also 
the central region of the bar where the Fe overabundance 
appears to be less prominent. Thus, observations seem to 
provide a consistent picture that the arc is where the rem- 
nant is interacting with dense ambient (possibly molecu- 
lar) clouds. On the other hand, the NW and SE parts of 
the SNR, where only the X-ray emitting ejecta are visible 
and there is no optical or IR emission, probably corre- 
sponds to low density regions where the forward shock in 
the ISM is still relatively fast and radiative shocks have 
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not formed. 

In fact, the projected location of SNR 0104-72.3 is 
at the eastern boun dary of the superbubble DEM SI 24 
(|Davies et al.l[T976h ( see Figure [3]). In Figure [3l^c), we 
show the Ha emission complex toward the direction of 
SNR 0104-72.3 and overlay locations of known SNRs in 
the SMC. It appears that SNR 0104-72.3 and SNR 0102- 
7219, a well-known core-collapse SNR, belong to the 
same star-forming region. The Type la origin inferred 
from the Fe-rich eject a and the association of the rem- 
nant with the star-forming region are not consistent with 
conventional scenarios of standard Type la, which involve 
the explosion of an old C/0 white dwarf in a binary. 
Recent studies of a large sample of Type la SNe sug- 
gested that progenitors of some Type la SNe ( "prompt" 
type) could have been relatively young massive stars 



with high metallicity (e.g. JScannapieco Bildstenl[2QQ5l : 
lAubourg et al.l l2008l ) . More recently, from the sys- 
tematic study of the SN Rs in the Magellanic Clouds, 
iMaoz fc BadenesI (|2010l ) suggested that some of the 
SNRs could be remnants of the prompt Type la SNe. 
While we cannot completely rule out the possibility of 
mere coincidence, the unusual nature of SNR 0104-72.3 
(i.e.. Type la characteristics of Fe-rich X-ray ejecta in 
the SNR and the SNR's spatial association with an en- 
vironment with a strong star-forming activity) suggests 
an intriguing possibility that SNR 0104-72.3 may be a 
candidate SNR of a prompt SN la. 

This work was supported in part by the Smithsonian 
Astrophysical Observatory under Chandra grant G08- 
9509 A. P.O.S. acknowledges partial support from NASA 
contract NAS8-03060. 
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